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ELECTROMAGNETIC  SCATTERING  BY 
CONDUCTING  BODIES  OF  REVOLUTION 

SOLUTION  USING 
SUB-  DOMAIN  AND  ENTIRE-DOMAIN 
BASIS  FUNCTIONS 


by 

J.  Joseph  and  R.  Mittra 
Electromagnetic  Communication  Laboratory 
University  of  Illinois 
Urbana,  Illinois 


Abstract 


In  this  work,  the  problem  of  electromagnetic  scattering  of  a  plane  wave  incident  on  a 
conducting  body  of  revolution  is  considered.  'The  body  is  assumed  to  be  situated  in  infinite 
homogeneous  space.  The  problem  is  solved  using  the  method  of  moments.  Use  of  two 
different  types  of  expansion  and  testing  functions,  namely,  sub-domain  type  and  entire- 
domain  type  is  considered.  Results  obtained  using  sub-domain  pulse  functions  arid  entire- 
domain  Gaussian  functions  are  presented.  The  relative  advantages  and  disadvantages  of 
each  type  of  basis  functions  is  discussed. 
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INTRODUCTION 


The  problem  of  elect  romagxiet  i  c  scattering  by  a 
conducting  body  of  revolution,  illuminated  by  a  plane  wave 
incident  from  some  arbitrary  direction  is  considered.  The 
geometry  of  the  problem  is  shown  in  Fig.l.  The  body  is 
generated  by  the  arc  C  revolving  around  the  z-axis . 

The  objective  is  to  compute  the  scattered  far  fields  and 
radar  cross-section  of  the  body  of  revolution.  In  this  paper, 
the  electric  field  integral  equation  which  enforces  the  total 
tangential  electric  fields  on  the  surface  of  the  body  to  be 
zero  is  used  to  solve  the  problem.  First,  equivalent  electric 
currents  are  postulated  on  the  surface  of  the  body.  The 
scattered  field  is  expresses  in  terms  of  these  equivalent 
currents.  The  electric  field  integral  equation  is  then  solved 
numerically  to  obtain  the  surface  currents.  The  surface 
currents  are  then  used  to  obtain  the  radar  cross-section  and 
far -fields , 

The  numerical  solution  of  the  integral  equation  is 
obtained  using  the  method  of  moments.  E'irst,  the  surface 
currents  are  expanded  in  a  set  of  appropriate  basis 
functions.  The  resulting  equation  is  then  tested  with  a  set 
o f  w e i g h t;  i n g  f u n c t i o n s  to  o b t ain  a  mat r :i .  x  e q u a t  i on  w i  t h  the 
coefficients  of  the  currents  as  unknowns . 

The  currents  are  expanded  in  a  Fourier  series  in  the 
az  imuthal  (<())  direction  and  tested  with  the  con  jugates  of  these 

expun. si  on  t  unct.  ions .  This  results  in  the  decoupling  of.  the 
•  •quat  i  ons  of:  various  azimuthal  mode:;.  The  advantage  of  this 
is  chat  one  can  solve  several  smaller  mat. rices  corresponding 


•cv" TV's. 


to  a  two  dimensional  problem  and  combine  these  solutions  to 
obtain  the  solution  for  the  body  of  revolution  problem. 

There  are  several  possible  choices  of  the 
expan s ion /test ing  functions  in  the  t-direct ion  (the  axial 
direction  -  see  Fig.  1)  for  the  current.  In  this  paper,  two 
classes  of  expansion  functions  namely,  sub-domain  and  entire- 
domain  functions  are  considered.  Two  programs,  PBOR  and  GBOR 
have  been  developed.  PBOR  uses  sub-domain  pulse  basis 
functions  and  testing  functions.  GBOR  uses  entire-domain 
Gaussian  basis  and  testing  functions.  The  advantages  and 
disadvantages  of  each  of  these  approaches  as  well  as  some 
numerical  examples  are  presented. 
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FORMULATION 


Consider  a  conducting  body  of  revolution  situated  in 
homogeneous  space  as  shown  in  Fig.  1.  The  body  is  illuminated 
by  a  plane  wave  incident  from  some  arbitrary  angle.  The 
total  electric  field  at  any  point  in  space  is  the  sum  of  the 
incident  field  and  the  field  scattered,  by  the  body.  The 
total  field  tangential  to  the  surface  of  the  conducting  body 
is  zero.  That  is. 


nxE 


total 


~  r _inc  ,  _scati  ^ 

=  nxlE  +E  J  =  0 


(1) 


on  the  surface  of  the  body.  In  (1),  Ei-nc  is  the  incident 
electric  field,  Escat  is  the  scattered  field,  and  n  is  the 
outward  pointing  normal  unit  vector  on  the  surface  of  the 
scatterer  (see  Fig.  1).  Using  equivalence  principle,  the 
scatterer  can  be  replaced  with  a  surface  electric  current  J 
radiating  in  homogeneous  space.  The  scattered  field  of  (1) 
can  then  be  expressed  in  terms  of  these  equivalent  currents 
as 

E'  ‘"(r)  =  - jC0A(r)-VO(r)  (2) 


with  the  potentials  of  (2)  given  by 

A(r)  =  —  f  j"  j(r')G(r,  r')dS'  (3  .) 

'UiJ  J 

C' 

o 

..  i  n  t'j 

4>  ~  —I  I  p(r')G(r,  r')dS'  (-1) 

A  71 1  J  J 


.vr- . 
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In  (3)  and  (4),  p.  and  e  are  the  parameters  of  the 

homogeneous  space  and  S  denotes  the  surface  of  the  scatterer. 
The  charge  density  p(r')  of  (4)  can  be  expressed  in  terms  of 

the  surface  current  J  as 

p(r')  =  i[V8,.j(r')J  (5) 


The  subscript  s  on  the  divergence  operator  denotes  surface 
divergence.  The  homogeneous  space  Greens  function  G  used  in 
(3)  and  (4)  is  defined  as 

—  jkR 

G(r,r')  =  — — —  (6) 


with  R  being  the  distance  between  an  observation  point 
r=(p,<t),z)  and  a  source  point  r  ’  =  (p  ’  ,  $  1  ,  z  ’  )  .  In  cylindrical 

coordinates  it  is  computed  as 


R  = 


r-r 


V  p  +p"~2pp  'cos  ($-<}>')+(  z- z  y 


(7) 


Now  we  can  combine  (1)  through  (7)  to  get.  an  integral 
equation  as  shown  below. 

nxEiru  -  nxj  I  j(r')c(r ,  r')dS'f-j— V  f  f  'V  J(r'))C(i ,  r')dS'j 

l  4  7tJ  J  47UO  J  J 


(8) 


Uur  -h  ji-ri  i  vr  now  is  to  no  1  ve  this  i  ntcqi  a  I  oquat  i  on  f  or 
!  he  oquiva  1  ent.  ;;ur  face  current,  J .  This  current,  can  be  used  to 
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Here  $  is  the  unit  vector  on  the  surface  of  the  body,  in 

the  cylindrical  coordinate  system  of  Fig.  1,  and  t  is  a  unit 
vector  on  the  surface  of  the  body  such  that  (n,<j),t)  form  a 

right  hand  triad.  The  incident  field  E^nc  of  (8)  also  can  be 
expressed  as  the  sum  of  its  <j>  and  t  components.  Thus  (8)  can 

be  written  in  component  operator  form  as 
"  E^nc 

T-,inc 

LE(j) 

The  Pij's  are  integro-dif f erential  operators  operating 
on  the  appropriate  components  of  current.  These  are  obtained 

by  computing  the  dot  product  of  (8)  with  appropriate  unit 
vectors  and  separating  out  the  part  and  J*-  part.  The 

explicit  expressions  for  these  operators  are  shown  below 


Pll  Pl2 
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(1U) 
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In  these  equations,  the  primed  coordinates  refer  to  the 
source  points  on  the  surface  of  the  scatterer  and  the 
unprimed  coordinates  refer  to  the  observation  points.  The 
angle  v  is  the  angle  tae  tangent  vector  t  makes  with  the  axis 

of  the  body.  This  angle  is  taken  to  be  positive  if  the 
tangent  vector  points  away  from  the  axis  of  the  body  and 
negative  if  it  points  towards  the  axis. The  surface  divergence 
in  surface  coordinates  given  by 


V-'-J  “  p'dt 


;t7r;  (  J &) 


pW 


has  been  used  to  obtain  (11) 


To  obtain  a  solution  for  the  unknown  currents  Jt  and  J<|) 
of  (11),  the  method  of  moments  is  used.  The  first  step  in 


the  solution  procedure  is  to  express  the  current  using  a  set 
of  suitable  basis  functions  as  shown  below. 


J(t,0)  -  0Jo(t,<j>)  +  tj,(t,4>) 
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•  ■*  +  y-xx-'-r.goe 
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Here  t  he  J  ’  are  the  unknown  coefficients  of  current  yet. 
e  it‘t  er  mined.  and  N^-  are  t.  he  number  of  basis 
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The  reason  for  solving  for  It11111  rather  than  will  be 
discussed  later. 

In  (12),  the  azimuthal  variation  of  the  currents  has 
been  expanded  in  a  Fourier  series.  The  basis  functions  used 
for  expansion  in  the  orthogonal  (t)  direction  are  denoted  by 
the  U's.  These  basis  functions  can  either  be  sub-domain  or 
entire-domain  basis  functions.  In  this  paper,  one  of  each 
kind  is  considered.  Sub-domain  pulse  basis  functions  with 
the  necessary  derivatives  approximated  using  finite 
differencing  and  entire-domain  Gaussian  basis  functions  are 
considered. 


Next  step  in  the  solution  process  is  to  test  the 
equations  with  a  set  of  suitable  weighting  functions  to 
obtain  a  matrix  equation.  Following  Galerkin's  method,  we 
test  with  the  conjugates  of  the  expansion  functions.  That 
is. 


rpPq 

=  T?e'iP*  =  u;(t)e'lp*', 

q=l 

r  '  t 

Nt 

(13a) 

rr.PS 

=  Tje-1^  =  Uj(t)e-Sp‘, 

q=J. 

r  r 

N» 

(13b) 

( 1  2  a )  is 

used  to  test  the  first 

row 

of 

(10)  and  (12b) 

is  used 

t  o  tet  t 

the  second  row. 

To 

obtain  the  expressions 

for 

t  he 

matrix  elements,  the 

Greens  f 

unction  G  appearing  in 

(11) 

i  s 

expanded  in  a 

Fourier 

se ries  as  shown  be  low . 
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where 

—  jkR 

Gm  =  cos(m^)  (14b) 

j- 71  R 

with  R  defined  in  (7)  with  ( <t> — <ji ’  )  replaced  with 

Testing  (10)  with  (13)  now  decouples  the  equations  for 
various  azimuthal  modes  m.  The  integrations  in  <|)  and  <|> '  can 

be  carried  out  analytically.  The  integrations  in  t  and  t' 
must  be  done  numerically.  Also,  the  Greens  function  Gm  is 

an  integral.  Thus  a  three  dimensional  numerical  integration 
is  required  to  obtain  the  matrix  elements  corresponding  to  p. 

The  generating  arc  C  in  Fig.l  over  which  the  t  and  t' 
integrations  are  to  be  performed  are  in  general  arbitrarily 
shaped.  For  computational  purposes,  this  curve  is 
approximated  as  linear  segments  connected  end  to  end  as  shown 
in  Fig.  2.  Then,  the  integrations  in  t  and  t*  must  be 
carried  out  over  each  of  these  segments.  To  reduce  the 
computational,  effort,  we  can  approximate  the  t  integral  by 
evaluating  the  integrand  only  on  one  point  per  linear  segment 
of  the  generating  arc.  This  is  equivalent  to  doing  point 
matching  when  the  testing  function  spans  only  one  linear 
segment,  such  as  is  the  case  when  pulse  basis  functions  are 
used.  Since  this  integrat  ion  is  carried  out  over  the 
incident  field  which  is  a  smoothly  varying  quantity,  this 
approximation  does  not  introduce  significant  error  if  the 
scatterer  is  modeled  with  sufficient  number  of  linear 
segments.  The  resulting  expressions  for  the  various  p’s  for 

azimuthal  mode  m  are  given  below.  The  first  summation  is 
over  the  field  segments  and  the  second  summation  is  over  the 
source  segments .  The  total,  number  of  linear  segments  is  F  and 
the  s'th  linear  segment  extends  from  t s_j  to  ts . 
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The  intf:grals  of  (15)  have  a  singularity  in  the 
integrand  when  the  source  point  coincides  with  the  field 
point.  For  a  discussion  of  the  handling  of  the  singular 
terms ,  see  [ 2 ]  . 

As  discussed  earlier,  the  basis  functions  U<j>  ana  U(-  can 

be  either  sub-domain  or  entire-domain  functions.  In  the  next 
sect  ic  i,  pulse  basis  functions  which  are  of  the  sub-domain 
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Sub-domain  pulse  basis  functions 


Fig.  3  shows  the  arrangement  of  pulse  basis  functions 
used  in  thi.s  paper.  P^'s  are  used  for  the  <j)  component  of 

current  and  '  s  are  used  for  the  t  component.  They  are 


defined  as  follows 


t,_n  ^  t  <,  t. 


uS(t)  =  P*(t)  = 


(l6a) 


otherwise 


S  t  <  t. 


U"(t)  =  P*(t)  = 


(16b) 


otherwise 


Observe  that  the  U0S  is  confined  to  the  s '  th  linear 

segment  .  Thus  there  are  a  total  of  F  of  these  basis 
functions.  That  is,  =  F .  On  the  other  hand,  we  have 

chosen  rhe  U^-s  to  straddle  over  two  linear  segments.  There 

are  a  total  of  F-l  of  these,  which  implies  that  =  F-l  . 


The  first  half  of  the  first  sub-section  and  the  last  half  of 


the  last  sub-section  is  assumed  to  have  zero  t-directed 


current.  This  is  not  a  correct  assumption  if  we  are 

considering  the  current  densities  for  a  scatterer  such  as  the 
one  in  Fig.  1.  However,  «re  can  solve  for  If-  =  2JTpJt,  rather 

than  Jt  itself.  In  this  case,  lt  will  be  close  to  zero  near 


the  points  where  the  generating  approaches  the  axis,  by 
virtue  of  p  approaching  zero.  This  is  the  reason  for 


expanding  the  currents  using  1^-  rather  than  using  Jt  in 
(12c)  . 


As  mentioned  previously,  the  testing  functions  used  are 
conjugates  of  the  expansion  functions.  Since  the  pulse  basis 
and  testing  functions  are  confined  to  a  small  region  of  the 
generating  arc  of  Fig.  2,  most  of  the  terms  of  the  double 
summation  of  (14)  will  drop  out.  The  resulting  expressions 
can  be  simplified  as  shown  in  (18)  .  In  obtaining  these 
expressions,  the  derivative  of  has  been  computed  using 

fini'e  differencing  as 


JLfn'Yt.p  - 

3tAUt^  0 . 5*(Atn+1+Atn) 


(17) 


As  mentioned  earlier,  the  t  integration  is  performed  by 
doing  one  function  evaluation  per  linear  segment.  This  turns 
out  to  be  one  function  evaluation  per  testing  function. 
These  function  evaluations  are  done  at  the  center  of  each 
testing  function.  The  expressions  for  the  elements  of  the 
moment  matrix  for  azimuthal  mode  m  using  pulse  basis 
functions  are  given  on  the  following  pages. 
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Entire-domain  Gau  ;ian  basis  functions 


Each  function  of  an  entire-domain  basis  set,  unlike  the 
sub-domain  basis  functions,  span  the  entire  structure.  The 
example  used  here  is  Gaussian  functions  as  shown  in  Fig.  4. 
The  n'th  basis  function  is  centered  at  Lln  (the  mean)  .  The 

"thickness"  of  the  basis  function  is  determined  by  its 
variance  <Tn.  The  mathematical  expressions  that  describe  the 

functions  of  Fig.  4  are  given  below. 

tL  <  t  <  t0 
0  <  t  <  tL 

t0  <  t  <  R 


(20) 


In  these  expressions,  the  quantity  R  is  the  length  of 
the  generating  arc.  t  g  and  tg  are  points  on  the  generating 

arc.  tjj  is  near  t-0  (the  beginning  of  the  generating  arc) 
and  tg  is  rear  t=R  (end  of  the  generating  arc)  .  The  basis 
functions  are  forced  to  go  to  zero  linearly  between  t j,  and  0, 
and  between  tg  and  R. 

The  basis  functions  used  for  both  Jt-  and  are  the  same 
and  are  given  in  (20)  .  However,  for  t  he  number  of  basis 
functions  used  is  two  more  than  that  used  for  Jt .  These  two 
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extra  basis  functions  are  also  Gaussian  basis  functions 
described  by  (20)  and  are  centered  at  t=0  and  t-P . 


The  derivatives  of  basis  and.  test~ng  functions  with 
respect  t  required  in  (15)  are  easily  computed  from  (20)  and 
are  given  below. 


tr,  <  t  <  trJ 


0  S  t  £  tL 


(  MR-OY 

e-°A - 


(-1) 

(R-tu) 


t0  £  t  <  R 


(21) 


Unlike  she  sub-domain  basis  set  where  the  number  of 
basis  functions  was  linked  to  the  number  of  linear  segments 
used  to  model  the  scatterer,  the  number  of  entire-domain 
basis  functions  used  have  no  relationship  to  the  number  of 
segments  on  the  generating  arc.  However,  since  the  t 
integrations  in  (15)  have  been  performed  by  evaluating  the 
integrand  at  one  point  per  linear  segment,  the  accuracy  of 
computations  are  dependent  on  the  number  of  linear  segments . 
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by  solving  a  smaller  matrix  than  would  be  required  if  sub¬ 
domain  basis  set  were  used. 


Incident  field 


An  arbitrarily  polarized  incident  plane  wave  field  Einc 
is  expressed  as  the  sum  of  two  linearly  polarized  plane  waves 
as 


/-.incAinc  .  —incli nc\  ~jkn  r 

(Ee  0  +  E0  0  )  e 


(22) 


(2  3  a) 
(2  3  b) 
(23c) 
(23c) 

I 

c. 

The  elements  of  the  forcing  function  vector  in  (10)  are 
^  found  by  finding  the  component  of  the  incident  field  of  (22) 

it  tangential  to  the  surface  of  the  scatterer  and  testing  with 

t h e  f  v n ct  ions  of  (13).  T he  r e s u  1 1  ing  ex p r e s s  i o n s  f  o  r  t;  h e 
S'  azimuthal  mode  m  are  shown  be  low .  I  n  these  expressions,  B[n 

refers  to  Bessel  fa net  ion  of  the  first  kind  of  order  m. 


C-'  whe  re 

.  nine  „  nine  line*  ,  Qinc  .  ,inc“  .  n^nc/' 

r'  0  =  cos0  cos<f»  x  +  cos0  sind,  y  --  sm0  z 

B  line  .inc-'  .  incA 

<J)  -  -“3J.n<p  x  +  cos(j)  y 

.  •  '  .  nine  line”  .  nine  ,  .ine*  ninc' 

■\  n  =  -sin0  cos<p  x  -  sm0  s .in<p  y  —  cos0  z 

r  =  pcos<p  x  +  psinp  y  +  zz 


♦I 
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*1 
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Summation  over  the  modes 

The  expressions  for  the  matrix  elements  and  forcing 
functions  given  in  (15)  and  (24)  are  for  a  single  azimuthal 
mode  m.  To  obtain  a  solution  to  the  scattering  problem,  one 
needs  to  sum  over  a  finite  number  of  modes  from  ~M  to  M.  The 
value  of  M  is  determined  by  several  factors  such  as  the  size 
of  the  scatterer  and  the  angle  of  incidence  of  the 
excitation.  Usually  the  fatter  (large  p)  a  scatterer,  the 

larger  the  number  of  modes  required.  When  the  excitation  is 
axially  incident,  the  only  modes  excited  are  1  and  -1.  The 
farther  the  excitation  is  away  from  the  axial  direction,  the 
more  the  number  of  modes  required  tc-  obtain  a  solution.  A 
rule  of  thumb  given  by  [3]  is 

M  ~~  kPmax-5in0inC+(5 

where  pmax  ■'■s  maximum  radial  dimension  of  the  scatterer. 

This  rule  is  valid  only  if  kpmaxsin9inc  is  greater  than  about 


The  expressions  for  the  matrix  elements  and  the  forcing 
functions  depend  on  the  mode  number  in.  This  dependence  is 
such  that  the  following  relationships  between  the  moment 
mat  r  ix  of  the  posit  ive  mode:.;  and  those  of  t  he  corresponding 
negat  i.  ve  modes  are  va  1  id  . 


I-  |V\ 


iv T  -iv; 


Liv.  iv\j  HV"  P."J 


Uiiniiai  iv,  i  lie  tore  i  mi  !  line!  ion:;  arc  rela'c,;  as  Iselb'ws  . 


v.V.v.-v  .Va 


s 


In  the  above,  the  superscripts  denote  the  component  of 
incident  field  and  the  subscripts  denote  local  coordinates  on 
the  scatterer.  This  relationship  between  the  forcing 
functions  and  matrix  elements  of  positive  and  negative  modes 
lead  to  the  following  relationship  between  the  coefficients 
of  equivalent  current. 


Vnc 
LJ0m  J 


Again,  the  superscripts  denote  to  the  component  of 
incident  field  that  is  producing  the  currents.  Because  of 
this  relationship  between  coefficients  of  current  for  the 
positive  and  negative  modes,  we  need  to  calculate  the 
currents  of  the  positive  modes. 


Computation  of  Radar  Cros3- section 

The  scattered  fields  can  be  computed  as  the  radiation 
field  field  from  the  equivalent  surface  currents .  Using 
reciprocity  [  i  },  the  expressions  for  this  field  can  be 
written  as 


E»u  = 


-j«H  e_jkR 


•J  ciS 


where  Er  is  the  field  due  to  a  test  source  located  at  the 
observation  point  with  its  amplitude  adjusted  to  produce  a 
unit  plane  wave  at  the  origin,  and  polarized  in  the  direction 
of  u.  Thus, 


For  the  moment  solution  where  current  is  given  as  a 
summation  of  the  product  of  the  basis  functions  Bn  and 

corresponding  coefficients  In,  (25)  can  be  vrritten  as 


E*u  = 


-jtOfl  e~jkR 
4k  R 


[r][i] 


(27  a) 


where 


(27b) 


the  integration  being  carried  out  over  the  surface  of  the 
scatterer.  In  the  above,  (27b)  has  exactly  the  same  form  as 

that:,  used  for  the  forcing  function  vector  discussed  earlier, 
the  only  difference  being  the  basis  functions  Bn  replacing 

the  testing  functions  Tq.  Since  we  are  using  Galerkin ' s 
procedure,  the  basis  and  testing  functions  are  conjugates  of 
one  another.  Hence  the  expressions  for  the  vector  [RJ  are 
obtained  from  the  expressions  for  the  forcing  function  given 
in  (24)  by  replacing  the  angles  0  nc  and  <|>  i n  G  with  the 

corresponding  observation  angles  and  changing  the  mode  number' 
m  to  "iit  to  account  for  the  conjugation. 


.  V  O  ■ 


Observe  that  in  (24)  we  have  broken  up  an  arbitrarily 
polarized  incident  plane  wave  into  two  linearly  polarized 
plane  waves.  In  terms  of  these  components  the  scattered 
field  components  can  be  written  as 


Eel  -jCOH  e~jkR  s60  s*l 


The  elements  of  the  scattering  matrix  are  obtained  .oy 
summing  the  contributions  due  to  the  various  azimuthal  modes 


-  I 


uv  being  00.  0<}>,  <J)0,  or  $<)).  The  terms  under  the  summation  are 


obtained  as  shown  below. 


=  K  J  [R*  mil 


K  mil. 


P",l 


PSl  PLJ  L  E0  m 


The  elements  of  the  various  R's  used  here  are  obtained 


from  the  expressions  foe.  the  E's  of  (24)  by  replacing  0inc 


a  i :  d  (J) 


w  i  L  h  t  h  e 


ipondi.nq  o b s e  r  v a  t.  i  o i i  a 1 1 g 


es  an-;: 


changing  the  mode  number  m  to  -m.  Once  t  he  r..  altered  field 


has  been  computed,  the  radar  cross-section  defined  as 


.■ .•*  oA'VA--'. ••V-V-Vh %’.•>  /-Vb'.v 


is  easily  computed.  Using  (28),  this  can  be  written  in  terms 
of  the  elements  of  the  scattering  matrix  as 


<3 


=  471  s 


This  can  easily  be  extended  ([4])  to  arbitrarily  polarized 
transmitters  and  receivers. 
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NUMERICAL  RESULTS 

In  this  section,  sample  results  obtained  using  the  sub- 
domain  basis  set  is  compared  with  those  obtained  using 
entire-domain  basis  set. 

First,  a  conducting  sphere  illuminated  by  a  plane  wave 
polarized  in  the  x  direction  and  travelling  in  the  -z 
direction  as  shown  in  Fig.  5a  is  considered.  The  radius  of 
the  sphere  is  1  wavelength.  The  pulse  basis  set  solution  is 
obtained  by  modeling  the  generating  arc  using  31  linear 
segments.  This  results  in  approximately  10  sub-sections  per 
wavelength.  (This  is  the  recommended  density  for  sub¬ 
sectional  basis  functions)  .  This  results  in  a  61  by  61 
matrix  equation.  For  the  entire  domain  basis  set,  the 
scatterer  is  modeled  using  the  same  number  of  linear  segments 

as  in  the  pulse  basis  set,  but  the  number  of  basis  functions 
used  is  11  for  Jj-  (and  13  for  J^)  .  This  is  approximately  4 

basis  functions  per  wavelength  and  is  the  recommended  density 
of  basis  functions  for  most  scatterers  when  using  entire 
domain  basis  functions. 

Since  the  excitation  in  this  case  is  axial,  the  only 
modes  excited  are  1  and  -1.  As  discussed  earlier,  only  the 
positive  numbered  nodes  are  computed,  since  the  solution  for 
the  correspond.!  ..g  negative  numbered  modes  can  be 
algebraically  derived  from  the  positive  mode  solutions. 

Fig  6  shows  the  comparison  of  bi. static  scattering 
cross-section  for  this  sphere  computed  using  the  two  methods. 
The  two  solutions  are  very  close  to  each,  other  .  Table  1 
compares  the  computation  times  required  for  each  method .  The 
time  used  by  the  entire- domain  program  is  aoout  twice  t hat 
required  for  the  sub- domain  program. 


Fig.  7  shows  the  comparison  of  bistatic  scattering 
cross-section  of  the  open-ended  cylinder  problem  shown  in 
Fig.  5b.  The  cylinder  is  4  wavelengths  long  and  the  matrix 
sizes  used  follow  the  rule  of  thumb  of  10  basis  functions  per 
wavelength  for  sub-domain  case  and  4  basis  functions  for  the 
entire  domain  case.  Again,  as  illustrated  by  the  numbers  in 
Table  1,  the  entire  domain  program  results  achieves  a  savings 
in  the  matrix  size,  but  at  the  expense  of  more  computer  time. 


DESCRIPTION  OF  PROGRAMS 


Tv-.'  programs,  referred  to  here  as  PBOR  and  GBOR  have 
been  developed  to  solve  the  conducting  body  of  revolution 
scattering  problem.  PBOR  uses  the  pulse  basis  set  and  GBOR 
uses  the  entire  domain  basis  set.  Both  program  require  the 
LINPACK  library  routines  CGECO, CGESL  and  CGEDI.  The  input 
data  is  read  from  a  file.  The  file  name  is  requested  by  the 
program  during  execution.  The  output  is  written  to  two  files 
whose  names  also  are  requested  by  the  program.  Notice  that 
the  input  data  required  by  the  two  programs  are  different. 

The  format  of  input  data  for  PBOR  is  shown  below. 

NFLDS ,  WVL 

THETA,  PHI,  ETHETA,  ETH-PHASE,  EPHI,  EPH-PHASE,  NOBS 

OBST, OBSP 


THETA,  PHI,  ETHETA,  ETH-PHASE,  EPHI,  EPH-PHASE,  NOBS 
OBST, OBSP 


RHO,  Z 

The  first  line  of  input  data  should  have  the  number  of 
different  excitations  (NFLDS)  and  the  wavelength  of  the 
incident  exci t at i on  ( s )  .  Following  this,  there  should  be 
NFLDS  lines  of  incident  field  data  with  each  line  having  the 
spherical  coordinate  angles  (in  degrees)  of  incidence  THETA 
and  PHI,  the  magnitude  (ETHETA)  and  phase  (ETH-PHASE)  in 
degrees,  of  t  he  t  heta  component  of  this  excitation,  the 

,  of  the  phi 


magnitude  (EPHI)  and  phase  (EPH-PHASE)  in  degrees 


component  of  the  excitation,  and  number  of  observation  points 
(NOBS)  at  which  the  far-fields  and  scattering  cross-section 
are  to  be  computed  for  this  excitation.  After  this,  there 
should  be  NOBS  lines  of  observation  point  data  with  each  line 
containing  the  observation  angles  theta  (OBST)  and  phi 
( OBSP )  .  Repeat  this  set  of  incident  field  data  and 
observation  point  data  as  many  times  as  there  are  incident 
excitations.  Finally,  there  should  be  as  many  lines  of  RHO 
and  Z  coordinate  points  as  required  to  describe  the 
scatterer.  These  points  are  connected  end  to  end  using 
straight  lines  to  obtain  an  approximation  to  the  generating 
arc  . 


The  format  of  the  input  data  set  for  GBOR  is  shown 
below . 


NFLDS ,  NOBS,  N,  WVL 

THETA,  PHI,  ETHETA,  ETH-PHASE,  EPHI,  EPH-PHASE 

OBST, OBSP 

RHO,  Z 


This  format  of  input  data  for  GBOR  is  somewhat  different 
from  that  of  PBOR.  The  number  N  indicates  the  number  of  basis 
functions  that  should  be  used.  It  should  be  such  that  there 
are  approximately  4  basis  functions  per  wavelength.  The 
number  of  observation  points  are  taken  to  be  the  same  for  all 
the  excitations.  Thus,  starting  at  the  second  line,  there 
should  be  NFLDS  lines  of  incident  field  data  and  following 
this  NOBS  lines  of  observation  point  data. 
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Table  1 .  Comparison  of  computation  times  and  matrix  sizes  for 

scattering  problems  using  entire  domain  basis  function 
and  sub-domain  basis  functions. 


